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Simulation  of  Life  Testing  Procedures  for 
Estimating  Long-Term  Degradation  and 
Lifetime  of  AlGaN/GaN  HEMTs 

Eric  R.  Heller 


Abstract — Finite  element  3-D  thermal  simulations  of  long-term 
degradation  in  AlGaN/GaN  HEMTs  for  high-power  applications 
are  reported  on,  in  which  temperature  evolves  over  time  as  the 
local  degradation  rate  varies  within  the  modeled  device  based  on 
the  local  temperature  of  the  degrading  region  (i.e.,  the  channel). 
Specifically,  hotter  regions  within  a  device  are  modeled  as  de¬ 
grading  faster  due  to  a  thermal  component  to  the  degradation 
rate  equation.  This  allows  self-consistent  simulation  of  life  testing, 
commonly  used  to  estimate  long-term  reliability  by  extrapolating 
failure  times  seen  at  elevated  channel  temperatures  to  a  lower 
“use”  temperature.  We  find  that  it  is  necessary  to  consider  the 
entire  distribution  of  temperatures  within  the  device  instead  of 
at  one  characteristic  location  to  get  the  most  accurate  estimates 
for  long-term  device  life.  The  effect  of  device  geometry,  assumed 
degradation  mode,  incorrect  thermal  resistance  data,  and  dissi¬ 
pated  power  level  on  this  lifetime  estimation  error  is  investigated. 
It  is  found  that  the  error  in  the  extrapolated  failure  time  is  greatly 
increased  when  both  the  thermal  resistance  is  in  error  and  the 
dissipated  power  of  the  life  test  does  not  match  the  expected  power 
during  operation,  compared  to  when  only  one  of  these  is  off. 

Index  Terms — Degradation,  FETs,  GaN,  GaN/AlGaN,  HEMTs, 
HFETs,  life  estimation,  life  testing,  MODFETs,  reliability,  simula¬ 
tion,  thermal  characterization,  thermal  resistance. 

1.  Introduction 

For  high-power-density  applications,  AlGaN/GaN 
HEMTs  have  been  identified  as  very  promising  de¬ 
vices  [1].  Because  of  the  high-power  capability  and  high- 
temperature  survivability  of  this  materials  system,  groups  have 
designed  devices  that  reach  very  high  channel  temperatures 
(Tch)  in  operation  [2]-[4].  Tch  can  easily  be  100  °C  or  more 
above  baseplate  temperature  and  can  vary  greatly  from  the 
center  of  a  central  finger  of  a  large  multifinger  device  to  the 
edge  [2]-[4]. 

A  common  method  of  estimating  the  long-term  operating 
lifetime  of  an  AlGaN/GaN  HEMT  for  power  radio-frequency 
(RE)  amplification  is  to  conduct  a  “life  test,”  where  several  (typ¬ 
ically  three)  populations  of  devices  are  run  at  different  baseplate 
temperatures  with  otherwise  similar  operating  conditions.  The 
devices  are  run  until  some  failure  criterion  is  reached,  which, 
for  a  dc  test,  is  usually  some  drop  such  as  10%  in  the  drain 


Manuscript  received  February  29,  2008;  revised  May  14,  2008.  Current 
version  published  September  24,  2008.  The  review  of  this  paper  was  arranged 
by  Editor  S.  Bandyopadhyay. 

The  author  is  with  the  Materials  and  Manufacturing  Directorate,  Air  Force 
Research  Laboratory,  Wright-Patterson  Air  Force  Base,  Dayton,  OH  45433 
USA  and  also  with  the  Semiconductor  Research  Center,  Wright  State  Univer¬ 
sity,  Dayton,  OH  45435  USA  (e-mail:  eric.heIIer@wpafb.af.miI). 

Digital  Object  Identiher  10. 1 109/TED. 2008.2003220 


X(Mm) 


Fig.  I.  Temperature  along  the  gates  (at  the  center  of  the  bottom  of  each  gate) 
for  a  representative  8  x  400  fim  device  with  IS-fim  spacing  between  gates 
running  at  Vds  =  10  V  and  Vqs  =  0  V,  and  temperature  at  the  baseplate  of 
200  °C.  Due  to  symmetry,  only  data  for  half  of  the  hngers  and  for  half  of  their 
length  are  shown. 

saturation  current  /dss  and,  for  RE  testing,  is  typically  a  certain 
drop  in  RE  power  or  operating  efficiency.  Higher  temperatures 
generally  lead  to  shorter  time  to  failure,  and  the  data  are  com¬ 
monly  fitted  to  a  thermal  Arrhenius  activation  energy  model 
for  the  failure  mode  responsible  for  the  degradation,  where  the 
temperature  required  in  the  model  is  that  of  the  failing  region 
[5].  This  Arrhenius  activation  energy  is  used  to  extrapolate 
the  measured  lifetimes  and  find  the  expected  lifetime  of  the 
device  at  the  baseplate  temperature  expected  during  operation, 
by  extrapolating  this  temperature  to  that  of  the  same  region  that 
was  known  to  fail. 

The  long-term  failure  of  a  good  device  is  commonly  at¬ 
tributed  to  degradation  over  the  entire  device  structure.  This 
may  lead  to  a  sudden  destructive  run-away  event  at  one  point, 
but  it  is  assumed  that  the  device  was  fairly  uniform  in  initial 
properties  at  the  start  of  operation  and  that  there  was  no  initially 
weak  point  that  drove  the  initial  degradation.  Of  course,  the 
device  is  not  at  a  single  uniform  temperature  during  the  high- 
power  operation  required  for  testing  these  devices.  Specifically, 
for  the  typical  large  device,  the  important  regions  of  the  device 
(i.e.,  channel,  gate,  source,  drain)  will  be  hottest  in  the  center  of 
the  device  and  cooler  at  the  edges,  as  shown  in  Eig.  1  [2]. 

Because  of  this  and  due  to  the  temperature-dependent  degra¬ 
dation  rate  assumed  in  the  Arrhenius  model,  the  local  degra¬ 
dation  rate  is  expected  to  vary  with  position  over  the  device. 
Traditional  failure  analysis  does  not  take  this  into  account  but 
uses  a  single  estimated  value  for  the  peak  temperature.  Device 
modeling  will  be  used  to  explore  how  the  data  extracted  from 
the  Arrhenius  model  changes  when  this  assumption  is  made 
compared  to  when  it  is  not.  As  time  progresses,  degradation 
is  seen  to  change  the  power  distribution  within  the  device,  and 


00I8-9383/$25.00  ©  2008  IEEE 

1 


Authorized  licensed  use  limited  to:  AFRL  Technical  Library.  Downloaded  on  June  15,  2009  at  13:01  from  IEEE  Xplore.  Restrictions  apply. 


HELLER:  SIMULATION  OL  LILE  TESTING  PROCEDURES 


2555 


both  this  and  the  resulting  evolution  of  temperature  with  time 
are  modeled.  In  addition,  temperature  varies  greatly  over  the 
source-gate-drain  (S-G-D)  region,  as  shown  in  an  earlier 
work  [2]  that  is  built  on  here  and  also  experimentally  by  Micro- 
Raman  in  [3].  Along  the  channel,  temperature  typically  peaks 
at  a  submicrometer-sized  region  at  the  drain  edge  of  the  gate 
where  the  device  is  pinched  off,  but  it  is  not  easy  to  get  this 
exact  temperature  [2]-[4].  It  is  expected  that  the  failing  region 
is  along  the  gates  where  the  temperature  is  highest,  and  so, 
this  temperature  is  used  for  the  failure  determination  in  this 
publication,  although  it  will  be  briefly  mentioned  how  data 
extracted  from  the  Arrhenius  model  change  when  the  region 
responsible  for  the  degradation  changes. 

Also  to  be  explored  is  the  change  in  the  thermal  conductivity 
with  temperature  of  some  of  the  materials  that  GaN  HEMTs  are 
composed  of,  specifically  SiC  and  GaN.  Of  course,  this  makes 
the  thermal  resistance  of  the  device  temperature  dependent,  and 
it  will  be  shown  through  device  modeling  that  this  temperature 
dependence  should  not  be  ignored  when  adjusting  the  baseplate 
temperature  during  a  life  test.  If  it  is,  it  can  lead  to  an  estimate 
of  the  activation  energy  that  is  considerably  higher  than  the  true 
value,  because  a  given  increase  in  the  baseplate  temperature  can 
lead  to  a  much  greater  increase  in  the  temperature  at  the  failure 
region.  Overestimating  the  activation  energy  for  failure  leads  to 
overly  optimistic  operating  lifetime  predictions. 

II.  Model 

To  model  these  effects,  a  representative  eight-finger,  400- fim 
gate  width,  and  25 -/im  gate  pitch  device  has  been  modeled  in 
the  finite  element  program  ANSYS  [6],  along  with  variations 
such  as  two  fingers  instead,  150-/im  gate  width,  35-/im  gate 
pitch,  etc.  The  GaN  was  1.2-/im  thick,  the  SiC  substrate  was 
350  yum  (or  175  yum  in  one  case),  and  the  SiC  simulation  domain 
is  2  mm  x  2  mm,  chosen  to  be  large  enough  so  that  further 
increases  do  not  change  the  model  results  significantly  [2],  [7]. 
All  boundary  conditions  are  adiabatic  except  for  the  bottom 
of  the  SiC  substrate,  which  is  fixed.  Because  of  fourfold  sym¬ 
metry,  only  one-fourth  of  the  device  is  modeled,  for  example, 
as  four  200-yum  fingers  with  adiabatic  boundary  conditions  at 
the  reflection  planes.  Because  of  degradation,  the  local  power 
dissipation  in  the  device  cannot  be  assumed  as  constant  along 
the  gate  width  but  is  assumed  to  vary  with  position  and  is 
loaded  into  an  initial  data  file.  Temperature-dependent  thermal 
conductivities  have  been  used  for  GaN  and  SiC  [8]  that  have 
been  previously  reported  in  the  literature  and  seen  in  our  prior 
modeling  efforts  [2]  to  reproduce  the  micro-Raman  results  of 
Kuball  et  al.  [3].  For  GaN,  this  prior  work  had  used  the  form 
in  [9],  but  for  this  work,  1.6  x  (T/293  K)“^  W/(cm  •  K)  was 
used  as  per  [3],  which  we  believe  may  better  reflect  the  tem¬ 
perature  dependence  of  the  GaN  thermal  conductivity  at  high 
temperatures.  Either  form  was  seen  to  give  the  same  qualitative 
conclusions.  Moreover,  it  is  known  that  the  interface  between 
GaN  and  SiC  will  have  a  low  thermal  conductivity  due  to  high 
defect  density,  also  known  as  a  thermal  boundary  resistance 
[10]-[12].  There  is  a  large  variation  in  the  reported  values 
of  this  resistance,  and  so,  for  simplicity,  this  effect  has  been 
neglected  here.  Fortunately,  it  is  not  the  absolute  temperature 


but  the  difference  in  temperatures  along  a  device  that  produces 
the  errors  investigated  in  this  model,  and  the  model  has  been 
found  fairly  insensitive  to  changes  of  this  nature. 

The  thermal  model  has  been  set  up  to  compute  and  store  the 
temperature  as  a  function  of  the  position  along  the  gate  width  of 
each  finger  for  this  arbitrary  power  distribution  data  file  that  was 
mentioned.  Instead  of  one  global  temperature  driving  degrada¬ 
tion,  it  is  assumed  that  there  is  one  point  within  a  plane  along 
the  S-G-D  cross  section  that  is  representative  of  the  position  of 
the  fail  region,  and  this  temperature  is  what  is  recorded  for  local 
degradation  rate  modeling.  The  power  distribution  data  (i.e., 
^DS  X  ^DS, local)  that  Were  used  by  ANSYS  and  the  temperature 
distribution  data  are  both  loaded  into  a  separate  utility,  and 
power  degradation  is  computed  by  a  user-defined  formula.  For 
all  of  the  simulations  that  were  run  for  this  work,  this  rate  was 
assumed  to  be  given  by 

P{x,  n,  t+At)=P{x,  n,  t)/(l+ At  (1) 

where  At  is  a  time  step  that  is  small  enough  that  the  degra¬ 
dation  rate  is  roughly  constant  during  its  duration  (i.e.,  the 
power  would  follow  a  decaying  exponential  with  time  if  T  was 
held  constant),  P  is  the  local  power  dissipation  at  one  place 
measured  along  the  gate  width  (which  extends  along  x)  for 
finger  n,  T  is  the  local  temperature  at  the  region  where  the 
failure  mode  is  assumed  to  take  place  (for  example,  the  center 
of  the  bottom  of  the  gate),  and  A  and  Ea  are  the  prefactor  and 
Arrhenius  activation  energy  for  the  degradation,  respectively. 
Data  for  this  time  step  are  archived,  time  is  incremented,  and 
the  power  distribution  table  is  recomputed  using  (1).  ANSYS 
can  now  be  used  to  compute  new  temperatures,  and  this  cycle 
is  repeated  to  the  end  point  time  where  the  degradation  has 
progressed  to  failure.  For  some  simulations  where  the  baseplate 
temperature  is  adjusted  during  the  run,  the  baseplate  input 
temperature  used  in  ANSYS  is  adjusted  once  per  time  step. 
Alternatively,  some  simulations  were  run  where  the  total  power 
was  held  constant  by  effectively  adjusting  the  gate  voltage;  this 
is  probably  more  common  than  is  baseplate  adjustment  for  dc 
life  tests  of  GaN  HEMT  devices.  In  this  case,  first,  (1)  is  used 
for  power  degradation,  and  then,  local  power  is  adjusted  to 
restore  the  total  power. 

Since  the  S-G-D  spacing  over  which  the  high  electric  fields 
exist  and  in  which  the  device  current  fiows  is  typically  ^  5  /im 
and  not  as  deep,  we  assume  for  our  modeling  effort  that  the 
physics  of  device  failure  (i.e.,  trap  generation  or  a  similar 
process)  operate  in  the  micrometer  or  smaller  size  scale,  and 
degradation  effects  at  one  place  in  the  device  do  not  directly 
affect  another  place  much  farther  away  on  this  finger  or  on 
another  finger.  For  reference,  the  gate  width  is  typically  ~  150- 
500  /im,  and  the  pitch  between  successive  gate  fingers  is 
typically  ^  25  /im.  Of  course,  there  are  indirect  effects  that  are 
included  in  the  modeling  effort.  Degradation  in  one  place  will 
affect  the  current  fiow  and  may  force  changes  on  the  gate  or 
drain  voltages  externally  applied  to  the  whole  device  depending 
on  the  testing  protocol  in  place.  Also,  local  power  dissipation 
creates  heat  that  diffuses  through  the  entire  device. 

The  ideal  power  distribution  will  have  slightly  less  power 
density  in  the  center  due  to  higher  device  temperatures  there. 
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Fig.  2.  Power  as  time  (hours)  progresses  for  an  8  x  400  initial  200  °C  baseplate  temperature  device  starting  with  an  initial  power  of  4.20  W/mm,  typical  of 
a  HEMT  running  at  Vds  =  10  V  and  Vqs  =  0  V  at  this  elevated  temperature.  For  clarity,  only  the  two  inner  hngers  and  the  two  outer  hngers  are  plotted.  In  (I), 
A  =  2  X  10®  and  Ea  =  1.5  eV,  and  the  temperature  used  for  Tfaii  is  at  the  center  of  the  bottom  of  the  gate. 


X  (pm)  X  (pm) 


Fig.  3.  Temperature  as  time  (hours)  progresses  for  the  same  device  and  simulated  test  as  in  Fig.  2,  with  the  temperature  at  the  center  of  the  bottom  of  each  gate 
plotted.  The  device  runs  from  —200  to  +200  fim.  Baseplate  temperature  was  adjusted  up  based  on  i^th  for  tho  center  of  the  inner  finger  at  f  =  0;  the  more  rapid 
degradation  in  the  device  center  causes  the  temperature  to  drop,  as  shown  there. 


We  expect  power  to  vary  as  ^1/T^{n}  (T  in  kelvins)  with 
n  ^  0.6  for  a  mostly  closed  channel  [2].  This  power  distribution 
was  modeled  and  found  to  not  be  a  significant  source  of  error, 
as  will  be  discussed.  Therefore,  for  simplicity,  the  initial  power 
density  has  been  assumed  constant  for  most  of  the  results 
reported. 

Last,  the  thermal  profile  within  a  device  can  vary  to  some 
degree  with  the  bias  used  to  achieve  a  given  power  dissipation. 
For  example,  a  fully  open  channel  device  has  heating  through 
the  entire  channel  and  also  from  contact  resistance,  while  heat¬ 
ing  for  a  pinched  off  channel  is  focused  on  the  pinch-off  region, 
even  when  the  total  power  is  conserved.  It  is  only  reasonable 
to  assume  that  degradation  can  affect  this  temperature  profile 
as  well,  as  will  small  changes  in  Vgs  to  adjust  the  total  power, 
both  of  which  will  then  affect  the  degradation  rate.  These  effects 
could  vary  with  the  physics  of  the  degradation  mechanism  and 
are  beyond  the  scope  of  this  paper. 


to  show  trends.  To  more  realistically  simulate  a  life  test,  the 
baseplate  temperature  has  been  continually  adjusted  up  to  try 
to  keep  the  channel  at  a  uniform  temperature,  as  may  be 
done  experimentally.  This  is  typically  done  by  measuring  or 
otherwise  determining  the  temperature  rise  over  the  baseplate 
temperature  for  the  desired  power  level  at  the  assumed  “failing” 
region  and  computing  a  “thermal  resistance”  i?th  by  dividing  by 
the  power.  In  practice,  this  can  only  be  done  approximately  be¬ 
cause,  as  was  shown,  the  device  has  a  considerable  temperature 
variation.  The  typical  way  to  approximate  this  experimentally 
is  to  compute  i?th  between  the  center  of  one  of  the  middle 
fingers  (the  hottest  place)  and  the  baseplate,  and  dividing  by 
the  total  power  dissipation  (Ptot)  of  the  device,  even  though 
this  is  not  the  same  as  the  local  power  dissipation.  The  failing 
region  (often  assumed  to  be  the  channel)  temperature  Tfaii  = 
Phase  +  Pth  X  Ptot  IS  then  held  constant  by  adjusting  Phase  as 
Ptot  decreases,  where 


III.  Results 

Fig.  1  shows  the  initial  temperature  distribution  along  the 
gates  of  the  simulated  device.  Figs.  2  and  3  show  the  power 
degradation  and  the  temperature  profiles  observed  when  a 
device  is  allowed  to  degrade  based  on  the  procedure  outlined 
earlier.  The  device  baseplate  temperature  is  200  °C  at  time 
{t)  =  0.  At  f  =  0,  the  full  initial  power  is  supplied  to  the 
device.  Temperatures  are  highest  in  the  center  of  the  centermost 
fingers,  and  so,  power  degradation  is  greatest  here.  Likewise, 
degradation  is  slowest  at  the  edges  of  the  device,  and  in  this 
example,  the  initial  degradation  at  the  far  edge  is  ^1.5%  of 
the  rate  in  the  center.  Degradation  has  been  allowed  to  progress 
for  a  ^19%  drop  in  total  power  (beyond  typical  failure  limits) 


Tfaii  =  Phase  +  Pth  X  Pfot-  (2) 

In  reality,  the  thermal  resistance  is  a  function  of  Phase,  Ptot, 
actual  device  bias,  local  power  dissipation,  etc.,  but  this  is 
often  neglected.  Fig.  3  shows  the  resulting  temperature  profile 
when  duplicating  this  procedure  with  the  model;  the  baseplate 
temperature  was  raised  from  200  °C  at  t  =  0  to  236.3  °C  at 
t  =  600  h  to  keep  Tfaii  constant  based  on  (2).  It  is  shown  that,  in 
reality,  the  temperature  varies  greatly,  even  in  the  center  of  the 
middle  finger,  and  in  this  case,  the  average  channel  temperature 
increased  by  17  °C. 

Fig.  4  shows  the  same  device  as  in  Fig.  2,  under  the  same 
initial  conditions.  In  this  case,  instead  of  adjusting  Tbase  for 
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Fig.  4.  Power  as  time  (hours)  progresses  for  a  device  similar  to  Fig.  2  and  for  the  same  initial  conditions,  except  that,  instead,  the  power  is  held  constant  by 
adjusting  Vqs.  and  baseplate  temperature  is  a  constant  200  °C. 


Fig.  5.  Temperature  as  time  (hours)  progresses  for  the  same  device  and  simulated  test  as  in  Fig.  4. 


Time  (hours) 

Fig.  6.  Degradation  of  average  device  power  shown  with  time  for  Ea  = 
1.5  eV  and  A  =  2  x  10®  with  200  °C  starting  baseplate  temperature.  Triangles 
show  the  degradation  based  on  the  modeled  temperature  profiles  of  Fig.  5, 
and  the  line  without  data  points  shows  the  rate  computed,  assuming  that  the 
entire  failure  region  is  at  Tf^ii.  The  thin  dash-dotted  line  shows  the  10%  power 
degradation  point. 

t  >  0  degradation,  Ptot  has  been  held  constant  by  adjusting 
Vgs  so  that  Tfaii  and  Tbase  both  remain  constant.  For  now,  it 
is  assumed  that  changes  in  Vgs  scale  all  portions  of  the  device 
up  or  down  in  power  by  the  same  percentage.  This  assumption 
will  be  discussed  in  more  detail  later.  Fig.  5  shows  the  resulting 
temperature  profile  for  Fig.  4.  The  final  temperature  device 
temperature  is  everywhere  cooler  than  before,  despite  the  same 
Pfaii  goal,  and  the  center  changes  far  more  than  before.  How¬ 
ever,  this  time,  the  average  channel  temperature  was  much  more 
constant  and  decreased  by  only  1  °C.  If  the  temperature  were 
monitored  by  some  means  that  extracts  an  averaged  temper¬ 
ature,  the  user  would  wrongly  think  the  temperature  did  not 
change  during  the  life  test. 

Fig.  6  shows  the  modeled  device  power  degradation  (line 
with  data  points)  based  on  the  temperature  profiles  in  Fig.  3 
versus  the  expected  degradation  when  it  is  assumed  that  Tfaii 


accurately  reflects  the  temperature  of  the  entire  failing  re¬ 
gion  as  is  commonly  done  (plain  line).  In  this  case,  Tfaii  = 
Tbase  +  Tth  X  Ptot  is  200  °C  +  0.04530  x  4200  =  390.3  °C. 
The  simple  one-temperature  model  yields  a  faster  degradation 
rate  because  Tfan  was  based  on  the  hottest  part  of  the  gate  finger 
(see  Fig.  3).  Of  course,  A  and  Ea  are  determined  based  on  the 
values  that  fit  the  observed  degradation  rate  and  not  the  other 
way  around.  The  point  of  Fig.  6  is  to  show  how  a  given  choice  of 
A  and  Ea  affects  the  simulated  degradation  and  that  A  and  Ea 
extracted  from  degradation  data  will  be  different  when  different 
channel  temperature  assumptions  are  made. 

As  mentioned,  the  ultimate  goal  of  the  life  test  is  to  find  the 
Arrhenius  activation  energy  Ea  and  prefactor  A  for  the  devices 
under  test,  and  this  cannot  be  done  exactly  when  the  temper¬ 
ature  of  the  failing  region  drifts  during  the  test.  To  determine 
the  extent  of  this  error  through  modeling,  the  simulated  life 
test  must  be  run  for  at  least  two  temperatures,  then  separately 
corrected  for  the  power  degradation  for  each  using  the  drop  in 
average  power  as  a  guide  (as  done  experimentally),  and  finally 
extract  E'^  and  A'  using  only  the  total  power  degradation  as 
must  be  done  experimentally.  In  all  cases,  10%  drop  in  total 
power  is  considered  a  failure.  These  quantities  are  primed  to 
differentiate  them  as  output  values,  separate  from  the  inputs 
used  in  (1).  In  this  case,  the  same  procedure  was  repeated 
at  Tbase  =  180  °C.  At  each  temperature,  the  “exact”  Pth  is 
extracted  from  the  thermal  model  at  the  device  finger  center 
at  t  =  0,  and  in  a  sense,  this  is  modeling  the  “ideal”  experiment 
without  measurement  errors  in  Pth-  Fig.  7(a)  and  (b)  shows 
the  results,  where  E'^  =  1.526  eV  and  A'  =  1.30  x  10^  are 
extracted  by  matching  the  time  for  10%  power  degradation  at 
this  temperature  and  simultaneously  for  Tbase  =  180  °C.  This 
contrasts  with  the  inputs  Ea  =  1.500  eV  and  A  =  2.00  x  10^. 
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Fig.  7.  (a)  Degradation  of  average  device  power  shown  with  time  after 

adjusting  Ea  and  A  for  the  curve  computed,  assuming  that  the  failure  region  is 
at  Tfaii.  Otherwise,  the  curves  are  the  same  as  in  Fig.  6.  (b)  is  the  same  as  (a) 
but  runs  with  a  lower  baseplate  temperature  of  180  °C  instead  of  200  °C.  The 
thin  dash-dotted  line  shows  the  10%  power  degradation  point. 


Two  temperatures  and  not  three  are  needed  to  extract  E'^  and  A' 
because  it  is  known  that  there  is  only  one  activation  energy  in 
the  modeled  degradation  rate.  A  few  test  cases  were  run  with  a 
third  temperature  to  verify  this. 

Simulations  where  the  attempt  to  hold  Tfaii  constant  was 
done  by  adjusting  gate  voltage  instead  of  increasing  Tbase 
are  reported  on  last,  because  of  a  complicating  factor.  This  is 
complicated  by  the  fact  that  the  proper  local  power  adjustment 
factor  depends  on  local  transconductance  (Gm).  which  may 
or  may  not  degrade.  Two  different  ways  for  the  device  to 
degrade  were  examined  so  that  it  would  be  known  if  this  was 
an  important  consideration.  In  the  first  scenario  (Scenario  1) 
investigated,  the  device  has  a  local  transconductance  {Gm, local) 
that  is  initially  uniform  and  also  does  not  drop  with  time. 
Degradation  in  this  case  is  instead  assumed  to  be  due  to  a 
locally  varying  threshold  voltage  shift.  In  this  case,  adjusting 
Vgs  simply  adds  the  same  value  everywhere  to  Ids, local  and 
therefore  also  adds  the  same  amount  of  power  to  the  local  power 
everywhere  in  the  device.  In  the  second  scenario  (Scenario  2), 
the  device  degrades  by  the  reduction  of  Gm, locals  but  local 
threshold  voltage  remains  fixed,  and  it  is  also  assumed  that 
the  entire  Vgs— G^, local  curve  for  all  Vgs  (at  least  near  the 
Vgs  of  interest)  has  degraded  by  the  same  percentage.  In  this 
case,  a  small  change  in  Vgs  will  everywhere  multiply  the 
^DS, local  and,  therefore,  the  local  power  by  the  same  factor. 
Effectively,  degraded  regions  have  less  power  recovery  when 
Vgs  is  adjusted  under  Scenario  2  than  under  Scenario  1.  Either 
way,  degradation  is  measured  by  monitoring  the  power  drop 
at  the  initial  gate  potential  and  is  equivalent  experimentally  to 
periodically  stopping  the  test  to  briefiy  return  the  gate  potential 


TABLE  I 

Extracted  E'^  and  A'  Estimates  When  Power  Is  Held 
Constant  by  Adjusting  Vqs  for  Scenario  I 


Rth  for  Tbase 
=180  °Crun 
(K-m/W) 

Rth  for  Tbase 
=200  °C  run 
(K-m/W) 

E'a 

(eV) 

A' 

Tfaii  for  Tbase 
=  180  °C 
run  (°C) 

Tfaii  for  Tbase 
=  200  °C 
run  (°C) 

0 

0 

1.250 

6.10  X  10^ 

180 

200 

0.015 

0.015 

1.613 

4.33  X  10“ 

243 

263 

0.03 

0.03 

2.023 

3.08  X  10’^ 

306 

326 

0.04245 

0.04530 

1.526 

1.17  X  10^ 

358.3 

390.3 

0.045 

0.045 

2.478 

2.19x1015 

369 

389 

TABLE  II 

Extracted  E'^  and  A'  Estimates  When  Power  Is  Held 
Constant  by  Adjusting  Vqs  for  Scenario  2 


Rth  for  Those 

=180  °Crun 
(K-m/W) 

Rth  for  Those 

=200  run 
(K-m/W) 

Eo 

(eV) 

A' 

TfaiAoY  Those 

=  180  °C 
run  (°C) 

Tfaii  for  Thase 

=  200  °C 
run  (°C) 

0 

0 

1.250 

5.78  X  10*1 

180 

200 

0.015 

0.015 

1.613 

4.11  X  10“ 

243 

263 

0.03 

0.03 

2.023 

2.92  X  1013 

306 

326 

0.04245 

0.04530 

1.527 

1.11  X  10^ 

358.3 

390.3 

0.045 

0.045 

2.479 

2.08  X  1015 

369 

389 

Time  (hours) 

Fig.  8.  Results  for  several  guesses  for  i?th  and  for  using  Vqs  adjustment  to 
keep  the  power  constant.  All  are  for  the  same  device  with  a  starting  baseplate 
temperature  of  200  °C.  The  thin  dash-dotted  line  shows  the  10%  power 
degradation  point. 


to  its  original  value  to  monitor  degradation  and  then  back  to 
resume  the  test  without  this  affecting  the  device. 

One  advantage  of  this  method  is  that  the  baseplate  tempera¬ 
ture  is  not  adjusted,  and  i?th  is  not  needed  during  the  life  test. 
i?th  is  still  needed  to  extract  E'^  and  A',  and  the  choice  of  Rth 
is  still  critical  for  correct  data.  Extracted  E'^  and  A'  are  shown 
in  Tables  I  and  II,  and  values  were  found  to  agree  by  better  than 
1  meV  and  10%,  respectively,  between  Scenarios  1  and  2  for 
the  same  i?th ;  these  differences  are  not  significant. 

Einally,  in  Eig.  8,  the  same  device  under  the  same  starting 
conditions  was  simulated  with  either  Tbase  adjustment  with 
an  assumed  Rth  or  Vgs  adjustment  to  keep  average  power 
constant.  In  practice,  Rth  is  not  an  exactly  known  parameter, 
and  it  is  common  in  a  practical  life  test  to  assume  a  fixed  value 
for  Rth  for  the  entire  test.  Eor  higher  Rth  values,  the  baseplate 
is  forced  to  greater  temperatures  to  keep  the  calculated  value  of 
Tfaii  constant,  and  this  is  why  the  degradation  rate  of  the  devices 
is  faster  for  these  runs.  Eor  Rth  =  0.045,  the  degradation  rate 
accelerates  slightly  with  time  despite  the  exponential  decay 
expected  from  (1)  because  Rth  is  too  high  and  the  modeled  tem¬ 
perature  profiles  (not  shown)  increase  with  time  everywhere, 
with  Tbase  =  218.9  °C  by  the  end  of  the  test.  Table  III  shows 
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TABLE  III 

Extracted  E'^  and  A'  Estimates  eor  Various  Assumptions  oe 
i^th-  In  All  Cases  Except  eor  i^th  =  0.  Baseplate  Temperature 
Is  Adjusted  in  an  Attempt  to  Hold  Tfaii  Constant.  For  All 
Scenarios,  the  Initial  “Correct”  Values  Used  in  the  Power 
Degradation  Computations  eor  Ea  and  A  Are  1.5  eV  and  2  x  10®. 
The  Fourth  Run  Does  so  Much  Better  Than  the  Rest  Because, 
IN  That  Case,  the  Known  Rth  From  the  Model  Was  Used 
IN  THE  Estimation  oe  Tfaii.  The  Region  Modeled  as 
Failing  Is  the  Center  oe  the  Bottom  oe  the  Gate 


Rth  for  Tbase 

=180  °C  run 
(Km/W) 

Rth  for  Tbase 

=200  °C  run 
(K-m/W) 

Ea 

(eV) 

A' 

Tfaii  for  Tbase 
=  180  °C 
run  (°C) 

Tfaii  for  Tbase 
=  200  °C 
run  (°C) 

0 

0 

1.264 

5.21  X  10‘^ 

180 

200 

0.015 

0.015 

1.607 

2.84  X  loll 

243 

263 

0.03 

0.03 

1.989 

1.47  X  1013 

306 

326 

0.04245 

0.04530 

1.526 

1.30  X  10^ 

358.3 

390.3 

0.045 

0.045 

2.413 

7.66  X  10i4 

369 

389 

TABLE  IV 

Extracted  E^  and  A'  Estimates  eor  Liee  Tests  oe  Dieeerent 
Devices.  The  “Baseline”  Device  Is  8  x  400  jum,  IS-jum  Gate 
Pitch,  Run  at  4.2  W/mm,  With  350-/im  SiC  Substrate,  and 
Liee-Tested  at  180  °C  and  200  °C,  With  Ea  and  A  Being  1.5  eV 
AND  2  X  10®,  Respectively.  All  Others  Are  Identical  to  the 
Baseline  Except  eor  the  Dieeerences  Listed  Here.  For  All, 
Power  Is  Held  Constant  by  Adjusting  Vqs.  Assuming  That 
Degradation  Proceeds  According  to  Scenario  2 


Device 

^a(eV) 

A' 

Z/aj/for  Tbase 

=  180  °C 
run  (°C) 

Tfaii  for  Tbase 

=  200  °C 
run  (°C) 

Baseline 

1.527 

1.11  X  10^ 

358.3 

390.3 

2.8W/mm 

1.518 

1.32  X  10^^ 

289.9 

317.1 

175  um  SiC  thk. 

1.525 

1.20  X  10^ 

330.3 

360.1 

2  X  400  fim 

1.505 

1.76  X  10^ 

269.1 

294.6 

2  X  150  |im 

1.504 

1.77  X  10^ 

257.9 

282.6 

Ea  =  2.0  eV 

2.035 

9.52  X  10^ 

358.3 

390.3 

35  fim  pitch 

1.521 

1.20  X  loll 

335.9 

366.2 

“Complex  initial 
power” 

1.527 

1.22  X  10^ 

354.9 

386.6 

the  results  for  E'^  and  A' \  it  is  shown  that  i?th  greatly  affects  the 
extracted  data.  Fig.  8  allows  for  a  comparison  to  Sozza  et  al.  [4], 
where  devices  were  stressed  under  constant  power  by  Vqs 
adjustment  (their  Fig.  6).  While  a  direct  comparison  cannot 
be  made  due  to  the  fact  that  this  is  a  different  device  under 
different  conditions  (this  paper  specifically  focuses  on  larger 
devices  than  theirs  with  greater  temperature  variations  within  a 
device),  it  can  be  seen  that  most  devices  have  a  degradation  rate 
that  slows  gradually  with  time  as  expected.  It  is  also  apparent 
that  there  can  still  be  a  large  device-to-device  variation  in  the 
degradation  rate  of  nominally  identical  AlGaN/GaN  HEMT 
devices,  which  makes  a  direct  comparison  difficult. 

Table  IV  shows  how  the  extracted  E'^  and  A'  vary  with 
variations  in  the  tested  device  geometry  and  conditions,  and 
can  be  used  as  a  guide  to  the  error  introduced  experimentally 
by  ignoring  the  temperature  variation  in  the  device  during 
a  life  test,  assuming.  In  other  words,  Ea  is  predicted  to  be 
overestimated  slightly,  and  A  is  underestimated  significantly 
compared  to  the  real  values.  As  expected,  variations  that  reduce 
the  heating  within  the  device,  through  either  reduced  power  or 
reduced  thermal  resistance,  will  make  the  extracted  parameters 
more  accurate.  This  is  particularly  true  for  reducing  the  ac¬ 
tive  area  of  the  device.  The  “complex  initial  power”  entry  of 
Table  IV  is  where  the  local  initial  power  was  set  as  P(x,  n,  t)  oc 


TABLE  V 

E'^  AND  A'  From  Ten  Liee  Tests  at  180  °C  and  200  °C  Are 
Extrapolated  to  Compute  Lieetime  eor  IOO  °C  Baseplate. 
The  200  °C  Baseplate  Device  Failed  at  378  h  and  the  180  °C 
Baseplate  Device  at  1465  h  at  4.20  W/mm,  and  7430  and  31  390  h 
AT  2.80  W/mm.  All  oe  the  Extrapolated  Lieetimes  (Last 
Column)  Should  Be  Compared  to  a  Simulation  oe  Degradation 
OE  THE  Same  Device  But  at  IOO  °C  Baseplate  Which 
Failed  at  1.50  x  10®  h  at  4.20  W/mm 


Rth  for  Tbase 
=180  °C  run 
(Km/W) 

Rth  for  Tbase 
=200  °C  run 
(Km/W) 

Power 
dissipation 
for  lifetest 
(W/mm) 

Ea 

(eV) 

A' 

Extrapolated 

tfaii  for  Tbase 

=100  °c 

0 

0 

4.20 

1.250 

5.78  X  10^ 

1.41  X  10^ 

0.015 

0.015 

4.20 

1.613 

4.11  X  lO^i 

1.14  X  10^ 

0.03 

0.03 

4.20 

2.023 

2.92  X  1013 

9.74  X  103 

0.04245 

0.04530 

4.20 

1.527 

1.11  X  10^^ 

1.51  X  10® 

0.045 

0.045 

4.20 

2.479 

2.08  X  1013 

8.62  X  103 

0 

0 

2.80 

1.331 

2.16  X  10^ 

4.71  X  10^ 

0.015 

0.015 

2.80 

1.584 

4.47  X  lOi® 

4.77  X  10^ 

0.03 

0.03 

2.80 

1.858 

9.23  X  loll 

6.68  X  103 

0.03925 

0.04182 

2.80 

1.518 

1.32  X  10^ 

1.05  X  10^ 

0.045 

0.045 

2.80 

2.154 

1.90X  10i3 

1.16  X  103 

T(x,  n,  to  account  for  the  expected  temperature  depen¬ 
dence  of  local  power  in  a  new  device  under  the  bias  conditions 
in  place,  as  mobility  and  saturation  velocity  drop  with  increas¬ 
ing  temperature  [2] .  This  relation  will  change  somewhat  under 
different  biases,  for  example,  with  the  exponent  increasing  to 
^0. 9-1.0  under  an  open  channel  [2],  [13].  As  an  example,  the 
initial  power  for  the  200  °C  baseplate  “baseline”  device  was 
4.085  W/mm  in  the  center  and  4.473  W/mm  at  the  farthest 
point,  but  with  the  average  value  still  4.2  W/mm.  Accounting 
for  this  effect,  the  180  °C  and  200  °C  devices  were  about  4° 
cooler  in  the  center  and  lasted  about  6%  longer  during  the  life 
test,  but  there  was  no  significant  effect  on  extracted  energy. 

The  same  modeling  process  used  for  the  life  tests  can  be 
used  to  simulate  the  device  life  under  operating  conditions.  In 
this  case,  the  same  degradation  mode  is  assumed  to  be  active, 
but  with  a  lower  baseplate  temperature.  It  is  also  assumed  that 
degradation  does  not  depend  directly  on  bias  conditions  or 
total  dissipated  power  but  only  indirectly  through  the  effect 
of  changing  temperature  in  (1).  Table  V  compares  lifetime 
calculated  using  E'^  and  A'  from  life  tests  at  180  °C  and  200  °C 
and  compared  to  modeled  lifetime  at  100  °C  baseplate,  where 
degradation  proceeds  according  to  Scenario  2.  It  is  shown  that 
the  extrapolated  lifetime  compares  fairly  well  to  the  modeled 
lifetime,  if  the  power  is  the  same  for  the  life  test  as  for  the 
simulated  operating  conditions.  This  is  true  even  when  E'^  and 
A'  are  off  because  the  wrong  Rth  was  used.  When  power  is 
different  for  the  life  test  versus  the  predicted  usage,  it  is  critical 
to  use  the  right  i?th-  The  observations  made  in  this  paragraph 
were  found  to  hold  true  even  if  fundamental  changes  are  made 
to  the  model  such  as  disabling  the  temperature  dependence  of 
the  thermal  conductivity  or  moving  the  region  where  Tfaii  is 
extracted  far  from  the  gate  (in  which  case  it  was  modeled  so 
that  both  the  failing  region  and  the  region  for  determining  Rth 
were  the  same  off-gate  region).  However,  if  a  different  region  is 
measured  to  get  Rth  through  error  or  limited  spatial  resolution 
of  the  measuring  technique  than  is  actually  failing,  then,  of 
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course,  E'^  will  be  off  as  discussed,  and  the  previous  statement 
does  not  apply. 

IV.  Conclusion 

It  can  be  seen  that  there  is  a  modest  error  introduced  in  the 
values  of  E'^  and  A'  extracted  from  a  life  test  from  neglecting 
the  spatial  and  variation  over  time  of  Rth  and  instead  using  the 
peak  value.  This  is  true  even  when  using  the  same  Rth  extracted 
from  the  thermal  model  for  the  simulated  life  test.  These  errors 
are  greatly  compounded  when  Rth  is  not  known  and  approx¬ 
imate  values  are  used  instead,  highlighting  the  importance  of 
getting  this  parameter  and  its  temperature  dependence  right  in 
a  life  test.  Of  course,  it  is  impossible  to  know  Rth  without 
knowing  exactly  where  the  fail  region  is,  and  even  then,  there 
may  be  considerable  uncertainty  in  this  parameter.  However, 
these  errors  can  be  minimized  by  careful  modeling  efforts  and 
by  careful  measurement  and,  in  general,  will  be  much  less  for 
smaller  devices. 

It  is  shown  that  the  error  in  E'^  and  A'  becomes  a  much 
greater  problem  for  lifetime  prediction  specifically  when  the 
dissipated  power  is  not  the  same  for  the  life  test  as  for  the 
usage  conditions.  Also,  it  is  seen  that  peak  and  average  fail 
region  temperature  evolve  differently  depending  on  how  the  life 
test  is  performed  (baseplate  versus  gate  voltage  adjustment),  so 
that  it  is  possible  for  this  average  temperature  to  stay  relatively 
constant  while  the  peak  temperature  changes  considerably;  an 
electrical  test  sensitive  to  this  average  temperature  would  be 
misleading. 

We  expect  the  same  issues  to  be  present  in  the  RF  device 
testing,  because  similar  temperature  profiles  as  discussed  here 
have  been  shown  experimentally  for  RF  drive  as  for  dc  drive 
[14],  and  temperature  is  again  an  accelerating  factor  in  degrada¬ 
tion  based  on  the  same  Arrhenius  model  (except  that,  typically, 
RF  output  power  replaces  Ids)-  However,  RF  life  testing  is 
commonly  done  by  driving  the  devices  to  a  set  power  saturation 
level  and  not  to  a  set  total  power  dissipation.  As  such,  there  is 
typically  a  different  amount  of  power  dissipation  during  the  life 
testing  than  there  is  in  the  device  during  actual  use,  and  so, 
based  on  our  findings,  errors  in  Rth  are  likely  to  be  much  more 
critical.  We  hope  to  extend  this  paper  to  allow  simulation  of  RF 
measurements  to  test  this  hypothesis. 
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